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Abstract 

The dynamics of soft (|p| ~ g 2 T) non-Abelian gauge fields at finite temperature is 
non-perturbative. The effective theory for the soft scale is determined by diagrams 
with external momenta po <g 2 T, |p| ~ g 2 T and loop momenta larger than g 2 T. We 
consider the polarization tensor beyond the hard thermal loop approximation, which 
accounts for loop momenta of order T. There are higher loop diagrams, involving also 
the scale gT, which are as important as the hard thermal loops. These higher loop 
contributions are characteristic for non-Abelian gauge theories and their calculation is 
simplified by using the hard thermal loop effective theory. Remarkably, the effective 
one-loop polarization tensor is found to be gauge fixing independent and transverse at 
leading order in the gauge coupling g. The transversality indicates that this approach 
leads to a gauge invariant effective theory. 
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1 Introduction 



Even at very high temperatures, when the running gauge coupling g is small, pertur- 
bation theory for non-Abelian gauge theories breaks down at the spatial momentum 
scale g 2 T [1], |2|. For static quantities, like the free energy or correlation lengths, one 
can integrate out the high momentum modes (p ^> g 2 T) in perturbation theory us- 
ing dimensional reduction |3|, [3J]. Then, at leading order, the non-perturbative physics 
associated with the scale g 2 T is described by a 3-dimensional Euclidean gauge theory, 
which can be easily treated in lattice simulations. 

Dynamical quantities, which are sensitive to the scale g 2 T, are more difficult to deal 
with. Since one has to consider non-zero real frequencies, one cannot use dimensional 
reduction. A prominent example for such a quantity is the rate for electroweak baryon 
number violation [], which, at leading order, is entirely due to gauge field modes with 
spatial momenta of order g 2 T ||. 

Fortunately, even for non-zero po, one can use perturbation theory to integrate out 
high momentum degrees of freedom to obtain an effective theory for the soft field 
modes. At leading logarithmic order, this effective theory is described by a Langevin 
equation 0. It determines the characteristic frequency of the soft field modes as 

p ~g*Tlog(l/g), (1.1) 

and it allows for a lattice calculation of the baryon number violation rate |§. 

The most direct approach to such an effective theory is to compute diagrams for soft 
external momenta, while the loop momenta are restricted to be larger than g 2 T. At 
one-loop order, the dominant contributions are the so-called hard thermal loops P|-|Tl|, 
which are saturated by loop momenta of order T. They lead to the Debye-screening of 
electric interactions on the length scale l/(gT). Therefore, only the magnetic degrees 
of freedom are non-perturbative. 

For a long time it was assumed, that the characteristic frequency of the soft mag- 
netic modes is of order g 2 T. It was realized in Ref. JT2J, that hard thermal loops 
drastically change this picture. They lead to a strong (Landau-) damping of the soft, 
non-perturbative dynamics. Considering the hard thermal loop resummed propagator, 
the characteristic frequency of the soft magnetic modes can be estimated as po ~ g 4 T 
0. 

In this paper, we investigate the polarization tensor for soft external momenta beyond 
the hard thermal loop approximation. There are higher loop contributions, which are 

2 For spatial vectors I use the notation k — |k|. 4-momentum vectors are denoted by — (k°, k), 
and the metric has the signature (H ) . 

3 More precisely, the Chern-Simons number diffusion rate, which is proportional to the baryon 
number violation rate close to thermal equilibrium H . 
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not suppressed relative to the hard thermal loops. They consist of self energy insertions 
on the hard loops and of ladder-type diagrams, and they can be economically calculated 
using the hard thermal loop effective theory. This paper describes the calculation of 
the effective one-loop polarization tensor, the leading logarithmic result was presented 
in [pl|. Summing the leading logarithmic contributions from all ra-loop diagrams by 
solving a Boltzmann equation for the soft field modes, one obtains the effective theory 
of Ref . |F| . Recently, this effective theory has also been obtained in , the Boltzmann 
equation was also obtained in |T3|-[16||. A recent detailed derivation of the Boltzmann 
equation can be found in |L7|. 

The paper is organized as follows: Sect. || briefly reviews some properties of hard 
thermal loops. Sect. [3] discusses higher loop contributions to the polarization tensor, 
which can be as large as the hard thermal loops. Sect. ^ contains the computation 
of the one-loop polarization tensor within the hard thermal loop effective theory in a 
covariant gauge. In Sect. [5] we discuss higher loop contributions and the relation of the 
diagrammatic approach to the Boltzmann equation. The results are summarized and 
discussed in Sect. |6|. Appendix A contains the expression for the sum over Matsubara 
frequencies which is used in the main text. Finally, it is shown in Appendix B, that 
the results of Sect. £| are also valid in Coulomb-like gauges. 



2 Hard thermal loops 



Hard thermal loops are one-loop contributions to n-point functions for external mo- 
menta small compared to T. They are saturated by loop momenta of order T. Fur- 
thermore, the dominant contribution is obtained when one momentum in the loop is 
on shell, Q 2 = 0. For the remaining propagators containing Q, one can use a large 
energy, or eikonal approximation 

1 1 1 (2.1) 



(Q + P) 2 2qvP ' 

where P is some linear combination of the external momenta, and t> M = Q^/q. The 
spatial part of v is the velocity of the hard particles in the rest frame of the plasma. 
Due to Eq. ( j2.1p , the angular integration over v and the integration over q factorize. 

The hard thermal loop 2-point function, which is the same in Abelian and non- 
Abelian theories, reads JTE| 



5n^(P) = mi 



+ PoJ- 



Air v ■ P 



(2.2) 



The q- integration is responsible for the factor m^, which is the leading order Debye 
mass squared. All fields coupling to the gauge fields contribute to the hard thermal 
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loops in the same functional form. In a SU(iV) gauge theory with Nf Dirac fermions 
and N s scalars in the fundamental representation one has = (1/3) (N + N s /2 + 
Nf/2)g 2 T 2 . Finally, the integral J dQ v is over the orientation of the unit vector v. 
Since ( |2.2|) is transverse, 

p^U^(P) = 0, (2.3) 

it can be written in terms of the 3-dimensionally (3-d) transverse and longitudinal 
projectors |1| Q 

7f(p) = VC°(p)=0, (2.4) 

vr(P) = ^-g^-rrip). (2.5) 



5W(P) = vr(p)sn t (p) + vr(P)snt(p). (2.6) 



such that 



For momenta Po^P^gT, the hard thermal loop ( [2.2| ) cannot be considered as a small 
correction to the tree level kinetic term, and it has to be resummed. In a covariant 
gauge with gauge fixing parameter £, the resummed propagator reads 



a^(p) = vr(p)\(P) + vr(P)MP) + (2.7) 

The 3-d transverse and longitudinal propagators 

A ' AP) = -p> + m tA py (2 ' 8) 

have poles at p^ = uj 2 e (p) > p 2 , and they have a cut for P 2 < 0. 

For static quantities, like the free energy or equal time correlation functions, it 
is convenient to work in the imaginary time formalism. Then po is an imaginary 
Matsubara frequency, and for p < gT the dominant contribution is given by the Po = 
modes. In this case the 3-d transverse fields are unaffected by hard thermal loops, 
<5n t (0, p) = 0, and the only effect of hard thermal loops is the Debye screening of the 
longitudinal fields due to SU^O, p) = m^. 

For dynamical quantities po has to be analytically continued towards the real axis. 
Different ways of approaching the real axis correspond to different time ordering pre- 
scriptions. For defmiteness, in this paper we will consider the case that 

Po = Re(po) + te, (2.9) 



Note that the projectors in Eqs. (2.4), (p.5[) have a sign different from the ones in Ref. pi 
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which gives the retarded propagator. 

When p Q is non-zero, also the magnetic, or 3-d transverse fields are screened ("dy- 
namical screening"). For p ~ g 2 T this has a dramatic effect [ 12| , since A t (P) is smaller 
than the free propagator by two powers of g when p ~ g 2 T. Therefore the transverse 
gauge fields perform only small fluctuations for this frequency scale. In order to obtain 
large fluctuations 5 , which are necessary for non-perturbative processes like electroweak 
baryon number violation, one has to consider the small frequency limit po <C p , in 
which case 

OT t (P)~-4^D- (po<p)- (2.10) 
4 p 

Then the hard thermal loop resummed propagator becomes 

A t (P)~l-^- (pb<p), (2.11) 

where 7 P = 4p 3 / (ftm 2 ^) . Since 7 P is of order g A T when p ~ g 2 T, Eq. ( |2.11| ) indicates 
that the characteristic frequency for the soft transverse modes is po ~ g A T The 
higher loop contributions discussed in this paper lead to a logarithmic correction to 
this estimate. 



3 Beyond hard thermal loops 

It will now become clear that the hard thermal loop approximation is not sufficient for 
obtaining the correct effective theory for the soft modes. Consider a hard thermal loop 
and imagine adding a self energy insertion on an internal line, where the additional 
loop momentum K is of order gT. This gives the diagram in Fig. |l|(a), in which the 
hard loop momentum Q is on shell, Q 2 = 0. Since the external momentum P is soft, 
the momentum Q + P is almost on shell. 

Compared to 511, the diagram in Fig. |l](a) contains two additional vertices, each of 
order gq ~ gT. There are two additional propagators containing the hard momentum 
Q, for which one can make the eikonal approximation (|2.1| ), and there is one hard 
thermal loop resummed propagator A(K) ~ l/(g 2 T 2 ). The loop integral over K 
comes with a Bose distribution function 

= e k /l_r (3 ' X) 

which can be approximated as n(k ) ~ T/k . Thus we can estimate 
5 For an instructive discussion, see Q . 
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(a) 



(b) 



Figure 1: Two-loop contribution to the polarization tensor. The external momentum 
P is soft (po, p<g 2 T) and the momentum Q is hard (q , q ~ T). The thick lines 
denote propagators with momenta of order T. The thin lines are hard thermal loop 
resummed gauge field propagators A(K) with k , k ~ gT. 



6U(P) x 



g 2 T 

v ■ P' 



(3.2) 



This shows that the diagram in Fig. [TJ(a) can be as large as 5U(P), when p Q and p are 
of order g 2 T or smaller. 

Now consider the diagram in Fig. [TJ(b) . Proceeding as above we estimate 



fl( lb )(P) ~ 5U(P) x (g 2 T 2 ^ ^ J d 4 k^- 



1 1 



TvKJ\Tv(K + P) 



\A(K). (3.3) 



On first sight, this diagram appears to be smaller than ( |3.2| ), since one eikonal prop- 
agator l/(y ■ P) got replaced by l/(v ■ K). However, performing the integral over ho, 
one obtains contributions from the poles of the propagators. For example, at the pole 
k = v ■ k, the second propagator in (|3.3| ) turns into l/(v ■ P). Therefore we can 
estimate 



II (lb) (P) 




5U(P) x 



v ■ P 



(3.4) 



which is as large as Eq. 

The above estimates for the diagrams in Fig. [I] are well known and have been dis- 
cussed extensively in the literature [0, |20], ^TJ. In order to deal with these large 
contributions, the right half of diagram (TJ(b) , together with the propagator A(K) has 
been interpreted as an insertion of a vertex correction. Then higher loop contributions 
were summed using a Schwinger-Dyson equation. For Abelian theories it was found 
that the large corrections cancel in the final answer. 
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00 

Figure 2: One- loop contributions to the polarization tensor in the hard thermal loop 
effective theory. The heavy dots are hard thermal loop vertices. Otherwise the notation 
is the same as in Fig. [I]. 

It should be emphasized that this view is not convenient, because it obscures the 
nature of the cancellation mechanism. Let us look at the diagrams from a different 
perspective ||13|| . Imagine first doing the integration over the hard momentum Q, with 



K kept fixed. This corresponds to integrating out hard momenta to obtain an effective 
4-point vertex for momenta small compared to T. Now we can make use of a well 
known fact: In an Abelian theory there are no hard thermal loop vertices for gauge 
fields only. In other words, the large contributions due to the diagrams in Fig. [I] cancel 
in exactly the same way as the diagrams which, by power counting, could give a hard 
thermal loop 4-point function. In particular, this means that the large contributions 
in an Abelian theory cancel even at a given loop order, i.e., to see this cancellation no 
resummation is necessary []. 

In contrast, in a no n- Abelian theory, there are hard thermal loop n-point functions 
for all n. Thus in non-Abelian theories this cancellation does not occur. 



4 One-loop polarization tensor within the hard ther- 
mal loop effective theory 

From the discussion in the previous section it should be clear that the calculation of the 
diagrams in Fig. |I| become simpler, if one uses the hard thermal loop effective theory. 
Then they correspond to a single diagram, which is shown in Fig. 0(a). 

The spatial loop momenta within the hard thermal loop effective theory should be 
smaller than A, where 

gT < A < T. (4.1) 

We also introduce a scale fi, which separates soft momenta from momenta of order gT, 
such that 

gy < // « #T. (4.2) 

6 The cancellation at fixed loop order has already been noticed in Ref. pM. 
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We consider only loop momenta larger than fi, which yields an effective theory for 
momenta smaller than [i. 

In addition to the propagator (|2.7|) , we will need the expressions for hard thermal 
loop 3- and 4-point functions. They are easily obtained from the hard thermal loop 
effective action in the form given in Refs. p3| . [] The n-point vertex can be written 
as 



n-2 2 I r '--v 

9 m B 



V-P n 



2tr(T an [T an - 1 ,[...,T 01 ]---]) 



v-P 1 v(P 1 + P 2 ) v(P l + --- + P n _ 2 ) 
+ permutations [(Pi, a x ), . . . , (P n _i, a n _i)] >, (4.3) 



where T a are the SU(iV) generators in the fundamental representation, normalized such 
that tr(T a T fe ) = (l/2)5 ab . 

We use the imaginary time formalism [|31|] . After the sum over Matsubara frequencies 
has been performed, one can analytically continue the external po towards the real axis. 



4.1 Diagram 2a 

With the expression for the hard thermal loop 4-point function taken from Eq. ( fO| 
one finds 



where 



k +Po 



k 



v-(K + P) v-K 



A pa (K), (4.4) 



fco=iu> n ' 



d 3 k 



f(K), 



(4.5) 



and u n = 2irnT are the Matsubara frequencies with integer n running from — oo to 
+oo. 

Performing the frequency sum can be quite tedious. One method, which turns out 
to be quite convenient in the present case, is described in Appendix A. In order make 
use of Eq. (|A.6| ), we rewrite the sum appearing in (|4.4| ) 



S 



K I 



kn 



kp+po 

v ■ (K + P) v-K 



A p(7 {K), 



(4.6) 



7 For other representations of the hard thermal loop effective action, see Rcfs. |2j-|27(]. A 

recent and particularly simple derivation of the effective action of Refs. J2J, [2?| can be found in p8f . 
For the relation of the different representations, see, e.g., 0, E^, |30||. 
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as 



s 



It 



k +Po 



k 



(K + P) v ■ K 



A pa {K + P) - A pa (K) 



Now we apply Eq. ( A.6 ) to obtain 
d 3 k 



(2tt) 3 
A^(po,k + p) 



(4.7) 



(4.8) 



A^(0,k) A p,7 (0,k + p 



A^(-p ,k) 



t> • P — v • k 

rfA; A^AT + \P) - A P °{K - \P) 



2-ni 



2P0) 



-p - v • k 

^0 + ho 



(K + \P) v (K - \P) 



(ko + §Po)(&o 

In this expression the analytic continuation of po away from Matsubara frequencies has 
already been performed. Furthermore, the high temperature approximation for the 
Bose distribution function has been used. 



Eq. (I~8) looks quite complicated. However, in the case, we are interested in, which 
is P ~ g 2 T, K ~ gT, it can be simplified considerably. Inside the integrand one can 
set P = except for the imaginary part of po- The error made with this approximation 
is of order g. Then, the first two terms in the curly bracket vanish, and one obtains 



„„„..,_ V" 11 ^ Tp °J 4vr (v-Py y(27r)*A%, V- -'J ■ \ r -K. 

Eq. (|4.9| ) is gauge fixing independent in covariant gauges. The terms proportional to 
the gauge fixing parameter £ contain factors v ■ K, which vanish due to 

1 



dQ v v^VuVpVa f d 4 k 1 



disc [A pa (K)) disc 



1 



• (4-9) 



disc 



v -K 



-i2n5(v ■ K). 



(4.10) 



Note that v ■ K — is the approximate on-shell condition for the loop momentum 
Q + P in Fig. |l|. In this case we have 



(vV t (k)v) = -(vV e (K)v) 



so that 



n^\P) * --mlNg 2 Tp,j 



dQ, 



1 _ 3i 

k 2 



{v-K = 0), 



(4.11) 



atx (v ■ py 



d A k 1 



disc ( A t (K) -A e (K)) disc 



1 



At this point one can see that is insensitive to the UV cutoff A. Due to Eq. 

( |4.10| ), only space-like momenta contribute to (|4.12|) . Then, the only discontinuity of 
the propagators A t £ is due to the discontinuity of the selfenergies, 



disc(A M (#)) = -Wi(K)\ 2 disc (SILtAK)) (A 2 < 0), 



(4.13) 
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Figure 3: Ladder-type diagram with two hard loop integrations corresponding to the 
effective diagram in Fig. 0(b). The notation is the same as in Fig. |l|. 



which falls off like 1/k 4 for k 3> mo- 

However, the part containing the transverse propagator it sensitive to the IR cutoff 
/i, because A t (i^) is unscreened for k <C k. To compute the /i-dependent piece of 
( |4.12| ), one can use the small frequency approximation ( |2.11| ) for A t (K) and neglect 
ko relative to k otherwise. Then the delta-function in Eq. ( |4 . 1 0| ) becomes <5(v • k), and 
one obtains 



Il^\P)^-^-mlNg*Tp J 



Ait {v-py 



log f — ] + finite 
V ^ J 



(4.14) 



where "finite" denotes terms which are IR finite, i.e., which are not sensitive to fi. The 
finite terms will depend on how the cutoff is implemented. 



4.2 Diagram 2b 



The result ( |1.14| ) is gauge fixing independent, but it is not transverse. There is, however, 
another one-loop diagram in the hard thermal loop effective theory, which contributes 
at the same order (Fig. H(b)). N In the original theory it corresponds to a 3- loop 
diagram with two hard loop momenta (Fig. |||). The vertices can be read off from Eq. 
~4[), and one obtains 



rf& pr {K)&. aX (K + P) 



dQ vl v^v lp v la [ dQ V2 v 2u v 2t v 2X 



Air 
k 



Vx ■ P J 4ir 
k +Po 



V\ ■ K v x -(K + P) 



v 2 -P 
k 



(4.15) 



k +Po 



v 2 -K v 2 -(K + P) 



Again, we use Eq. ( |A.6| ) to perform the sum over Matsubara frequencies, and the terms 
f(0,po), f(po,Po) in ( |A.6p do not contribute at leading order in g. We neglect P relative 

8 The diagrams in Fig. ^ have been evaluated previously on the plasmon "mass shell" of the prop- 
agators ( |2.8|) to compute the gluon damping rate |Q and corrections to the longitudinal plasmon 
frequency J33|| . Our result does not apply to this case since we consider po <C gT. 



9 



to K except for the imaginary part of po (see below), which gives 



rev) 



--m D Ng Tp 

d A k 



dtt Vl Vlf.VlpVlg f dn v2 V 2v V 2t V 2 \ 
47T Vl • P J 4:71 



V 2 -P 
1 



(4.16) 



J——A^(k - ie, k)A° x (k + ie, k)disc (^-^) disc 



1 



v 2 -K 



Due to the delta-functions 8(vi ■ K), only space-like K contribute in Eq. ( 4.16Q . Thus, 
there is no contribution from the plasmon poles of the propagators at k\ = uj 2 e (k) > k 2 , 
which means that large denominators as in Eq. ( |2.1| ) cannot arise here. Therefore, it 
was consistent to neglect P in the propagators. 

Now we insert the expression for the propagator ( |2.7| ). Again, the terms proportional 
to the gauge fixing parameter £ drop out due to the delta functions 5{v ■ K). Thus, 
we find that also the leading order contribution from diagram 2(b) is gauge fixing 
independent. The remaining terms give 

, 'IP) ~ --in! \ <r'ln,, I 1- li— / 

[1U 



rt£\P) 



2™ B N g T Po j—— 



d 4 k 



V2u 



Atx v 2 ■ P 



(4.17) 



(2tt)< 



{v{P t v 2 ) A t {K) + {v{P t v 2 ) A e {K) 



disc 



vi ■ K 



disc 



1 



v 2 -K 



It is now easy to see that the i^-integral is insensitive to A: for k 3> gT the propagators 
fall off like 1/k 2 . Note that this integral also appears in the collision term in the 



Boltzmann equation of Refs. [fL4], [17| . There, the integrand is given by the square of 
an on-shell matrix element. 

Again, there is an infrared sensitive piece in the i^-integral, which is due to the 
term containing |A t (i^)| 2 . To compute this piece, one can use the small frequency 
approximation ( |2.11| ) for A t (K), and one can neglect ko otherwise. Using 

2 ( Vl ■ v 2 ) 2 



dn k (viV t v 2 y5(vi-k)5(v 2 -k) 



k 2 



1 



vVl • v 2 



ei-K' 



(4.18) 



one finds 



ngb) (P ) 



—m D Ng Tp 



log 



'gT 



dttvi Vi^ fdVt V2 v 2 
4ir vi- P 
(vi • v 2 ) 2 



1 - (Vi • v 2 ) 



47T V 2 ■ P 

+ finite . 



(4.19) 



4.3 Diagrams containing tree level vertices 

So far we have considered diagrams, which contain only hard thermal loop vertices. 
Now we will see that diagrams containing tree level vertices are smaller by one power 
of g. 
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Figure 4: Sub-leading contribution to the one-loop polarization tensor with one hard 
thermal loop vertex and one tree level vertex. The notation is the same as in Fig. |2|. 

The diagrams containing only tree level vertices have already been discussed in 
Ref. J7j. They are quite similar to hard thermal loops since the loop momentum is 
large compared to the external momentum. The size of hard thermal loops is deter- 
mined by integrals like 



dkkn(k) ~ T 2 



(4.20) 



In our case, k is smaller than A, which means that one can approximate n(k) 
Thus, instead of ( 4.20|) , we have 



T/k. 



dkk 



k 



(4.21; 



The integration region k ~ gT gives a contribution which is smaller than the hard 
thermal loop by one power of g. From k near the cutoff A, one obtains a contribution 
which cancels the A-dependence of the integral ( f4.20|) 

Now consider the diagram with one hard thermal loop and one tree level vertex in 
Fig. [|, which gives 



ng(P) 



i 2 fdQ v v fl v p v a 



^ A^(K)A ffX (K + P)T uXt (-P, K + P, -K) 



k fcp +p 

v ■ K v(K + P) 



(4.22) 



where T is the tree level vertex. After performing the frequency sum, the square bracket 
gives a factor 5(v ■ K) in the P — > limit. Thus, there are only contributions from 
space-like K and one can also take the limit P — >• in the propagators. Since fl4.22|) 
contains only a single factor 1/v ■ P rather than two, it is suppressed relative to ( f4.14|) , 
( |4.19| ) by a factor g. 

9 The difficulties with this cancellation which were encountered in Ref. J34| arc irrelevant as long 
as one does perturbation theory. 
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4.4 The logarithmic approximation 

Keeping only the terms in Eqs. ( |4.14j ), ( |4.19| ), which depend logarithmically on fi, the 
polarization tensor becomes 

n^^W.sjfj./^i/^^.K,,. (4,3, 
Here, the function 

J( Vl) v 2 ) = -^(vi - v 2 ) + - 2 (Vl ' V2)2 (4.24) 

7T 2 _ ( Vl . V2 )2 

is the same as the one which appears in the noise correlator and in the collision term in 
the Boltzmann equation of Ref. (cf. Sec. |5|). Furthermore, 5^ s2 ^ is the delta-function 
on the two dimensional unit sphere, 

|^ v J(v 1 )5 (52) (v-v 1 ) = /(v). (4.25) 



One can easily verify that 



J dQ Vl 7( Vl ,v 2 ) = 0. (4.26) 



This implies that the polarization tensor ([4.23|) is transverse with respect to the external 
4-momentum, 

p^ u A \P) = 0. (4.27) 

This condition is necessary for the effective theory for the soft fields to be gauge in- 
variant. 



4.5 Beyond the logarithmic approximation 

In the previous section we found that the leading logarithmic result for the polariza- 
tion tensor is transverse. Now we will see that this result holds beyond the leading 
logarithmic approximation, i.e., at leading order in g. 

In order to define IL^ beyond the logarithmic approximation, one has to specify 
the way how the loop momenta are cut off in the infrared. A convenient method is 
dimensional regularization. Then, the above calculation of 11^ has to be repeated in 
n = 4 — 2e instead of 4 dimensions, which is straightforward. In the hard thermal loop 
selfenergies and vertices, the angular integration J dQ(- ■ •) become n — 2 dimensional, 
and in Eqs. ( |2.2|) , (fO|) one has to replace 




(4.28) 
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where fi = / dVt. Furthermore, one has to introduce a scale p, in order to keep the 
coupling constant dimensionless. 

Then, the polarization tensor can be written as 



TT / ry\ ■ 2 ]\T 2-2erp f «^Vl f <^v 2 V 2u f, \ 



tt Vl - P J tt V 2 ■ P 



where 
-T(vi, v 2 ) 



1 /• rf n fc 



2 7 (27r) n 
'^"- 2 )( Vl -v 2 )l 

Kq 



(4.29) 



(4.30) 



/c 2 



disc 



A t (K) - A e (K) 



{v{P t v 2 ) A t {K) + (v{P e v 2 ) A e (K) 



disc 



vi ■ K 



disc 
disc 



vi-K 
1 



v 2 -K 



We will now see, that /(vi, v 2 ) also satisfies ( |4.26| ), which then implies that IT in Eq. 
Q4.29Q is transverse. We have 



dQ 



Vl 



n 



J(vi,v 2 ) = 

i f 1 - f) disc 



1 r d n k 



2 J (2 



-disc 



7T 



v 2 -K 



A t (K)-A £ {K) 



(4.31) 



Vl 



n 



[viP t v 2 ) A t (K) + (viV e v 2 ) A e (K) 



disc 



■ K 



For the first term in the curly bracket one can use Eq. (|4.13| ) to write the discontinuities 
of the propagators in terms of the discontinuities of the corresponding selfenergies. The 
latter can be read off from (the n-dimensional version of) Eq. Q2.2|) making use of the 
projectors (|2.4| ), fl2.5|) . Taking into account that, in n dimensions, V[ l = n — 2, one 
finds 



disc 



A t {K)-At{K) 



-mlkn 11 - M 



dtt 



k 2 



n-2 



\A t (K)\ 2 + \A e (K)f 



Vl 



n 



disc 



v x ■ K 



(4.32) 



which contains precisely the same Q Vl -integral as the second term in the curly bracket 
of Eq. fl£3Ip . 

In the second term in the curly bracket of Eq. ( [4.31] ) we can perform the angular 
integration over vi. Due to Vi ■ K = v 2 ■ K = 0, the velocity vectors can be written as 



I- 2 



.2 ^ 



(i = l,2), 



(4.33) 
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where (pi are unit vectors orthogonal to k = h/k. Then we have 



{viVtV 2 ) 



1 - 



■1 + 



k\ 
k 2 

k 2 ' 



(4.34) 
(4.35) 



Now one can perform the y3i-integration. The transverse-longitudinal interference term 
in Eq. ( |4.31| ) vanishes due to 



J d0i0i -02 = 0. 



Furthermore, we have 



dipl {0! ■ (pi 



n-2 



(4.36) 



(4.37) 



Consequently, for t> 2 • K = 0, one can write 



da 



n 



[vtV t v 2 ) A t (AT) + (viPtVz) At(K) 



disc 



v\ ■ K 



1 k 2 



Vl 



n 



n-2 



\A t (K)\ 2 + \A e (K)f 



disc 



vi ■ K 



(4.38) 



Inserting Eqs. ( f4.32p , ( OS ) into ([4.31 ), one finds that the two terms in the curly 
bracket cancel, so that 



dft Vl /(vi,v 2 ) = 0. (4.39) 
This implies that II in Eq. (|4.29 ) is transverse with respect to the external n- momentum, 

p^ u (P) = 0. (4.40) 



5 Higher loops 

In the previous section we have seen that there are one loop diagrams which are as 
important for the soft field modes as the hard thermal loop 2-point function. The 
natural question arises whether there are higher loop diagrams which contribute at the 
same order and which have to be included in the effective theory for the soft fields as 
well. 

The computation of higher loop diagrams along the lines of Sec. 4 can be expected 
to be quite tedious. The number of terms in the hard thermal loop n-point functions 
( fO| ) grows rapidly with n. Moreover, performing the sum over Matsubara frequencies 
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in multi-loop diagrams is probably very difficult (even for the one-loop diagram in 
Fig. 2b it was not particularly easy). A key simplification made in this paper was the 
high temperature approximation for the Bose distribution function, which corresponds 
to the classical field limit. Unfortunately, this simplification could be made only after 
performing the Matsubara sum. Therefore it is more convenient to use an alternative 
formulation of the hard thermal loop effective theory, which incorporates the classical 
field approximation right from the start. 

Fortunately, such a formulation exists [22, 2^J. It is the non-Abelian generalization 
of the linearized Vlasov equations for an electro-magnetic plasma (see e.g. |35| ). In 
addition to the classical gauge fields, these equations contain fields W a (x, v) describing 
the fluctuations of the phase space density of particles with momenta q of order T and 
q/q = v. The W a transform under the adjoint representation of the gauge group. In 
this formulation it is also possible to integrate out momenta of order gT perturbatively 
|7], 36 1 . Again, one encounters terms which are logarithmically sensitive to the se- 



paration scale /i. Keeping only the leading logarithmic terms one obtains an effective 
theory which is described by the Boltzmann equation 0, 

v-DW{x,v) = v-E(x)+£(x,v)- f^^C(v,V)W(x,v r ), (5.1) 

J 4:71 

together with the the Maxwell- Yang-Mills equation 

D^ v {x) =ml[ ^v v W(x, v). (5.2) 

£ is a Gaussian white noise which is specified by its 2-point function 

IT 

(C(*,v)£V,v')> = ^ 4 (x-xWC(v,v'). (5.3) 



m 



The integral kernel of the collision term reads 

C(v,v') = —Ng 2 T log I(v,V), (5.4) 

where /(v, v') is given by Eq. (|4.26 ). Correlation functions for soft external momenta 
are obtained by solving ( |5.1| ), ( |5.2| ) and performing the thermal average over initial 
conditions together with the noise average. 

We will now make contact between the Boltzmann equation ( |5.1| ) and the one-loop 
calculation in Sec. [|. This will demonstrate that there are indeed higher loop contri- 
butions to the effective theory for soft field modes, which are summed by using the 
Boltzmann equation. 

Loosely speaking, the lhs of Eq. ( |5.2|) is the first functional derivative of the effec- 
tive action which contains the effect of the field modes with momenta larger than fi. 
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Now we consider the analogous quantity for the diagrammatic approach. We denote 
the generating functional of the the "tree level" contribution (|2.2|) plus the one loop 
contribution (|4.23|) by T. The first (functional) derivative of T is 



ST[A] 



5A V (-P) 

Obviously, this can be written in the form 



m» v {p) + Yi% A \p)] a»(p) 



(5.5) 



ft[A] 
5A"{-P) 



ITLt 



v v u W(P,v) 



Alt 



where 



W = W + W 2 . 



(5.6) 



(5.7) 



Here the subscripts 0, 2 count the powers of g which are not contained in m^. Wo, 
corresponding to SH, reads 



Wo = -A + v ■ A, 



v ■ P 



(5.8) 



and W 2 , representing H^ U A \ is given by 



W 2 (v) 



dQ„' 1 



V ■ P J 47T V' ■ P 



C(v, v')v' ■ A. 



Eq. ( |5.9| ) can also be written as 



W 2 (y) = -% 



1 f dQ v i 



v ■ P J Air 



c(vy)w (v'), 



(5.9) 



(5.10) 



because due to Eq. ( |4.26|) the term —A in W drops out after integrating over v'. W 
satisfies the usual linear Vlasov equation 



- iv ■ P Wo = iv l + p°A l = (v • E) linear . 



For W 2 we have 



-iv-P W 2 {v) = - / ^C(v, v') W (V). 

Air 



Adding (|5.11|) and (|5.12j ) we obtain 

-iv-PW(v) = (v-E) linear 



Air 



C(v,v')W (V) 



(5.11) 



(5.12) 



(5.13) 



Now we see that W satisfies the linearized Boltzmann equation (without noise) up 
to terms which are higher order in the loop expansion parameter g 2 . Thus, from the 
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validity of Eq. ( |5.1| ) one can infer that there are indeed higher loop contributions to 
the leading order effective theory for the soft field modes. Furthermore, the fact that 



the lhs of ( |5.1| ) contains a non-linear term implies that there are also vertex functions 
which have to be included in the effective theory for the soft fields. 

It is interesting to note that by writing the one loop result computed in this paper 
as in Eq. ( |5.13|) , one can already "guess" the form of the Boltzmann equation. All 
one has to do to get from Eq. ( |5.13| ) to ( |5.1| ) is to replace Wo by W and to include 
non-linear terms to make the equation gauge covariant. However, in this way one does 
not obtain a noise term. A more quantitative discussion which accounts for the role of 
the noise will be deferred to a separate publication. 



6 Summary and Discussion 

We have seen that in non-Abelian gauge theories the hard thermal loop approximation 
is not sufficient for obtaining the effective theory for soft (p ~ g 2 T) momentum fields. 
When po < g 2 T, the polarization tensor for soft external momenta receives higher loop 
contributions from loop momenta larger than g 2 T, which are as large as the hard 
thermal loops. 

These large higher loop contributions are due to ladder type diagrams and diagrams 
with selfenergy insertions on propagators carrying hard momenta. Their calculation is 
simplified if it is performed in two steps. First, one integrates out the hard momenta, 
which gives the well known hard thermal loop effective theory for momenta p T. 
Then, the polarization tensor is calculated within this effective theory by integrating 
over momenta of order gT. When the calculation is organized in this way, it is easy to 
see why these large contributions cancel in an Abelian theory. The large contributions 
discussed in this paper contain only hard thermal loop gauge field vertices, which are 
absent in the Abelian case. 

The calculation of the one-loop polarization tensor in the hard thermal loop effective 
theory requires an explicit IR cutoff separating spatial momenta of order g 2 T from 
momenta of order gT. We have obtained an expression for the polarization tensor valid 
to leading order in g, which is entirely due to diagrams containing only hard thermal 
loop vertices. They correspond to a sum of 2- and 3-loop diagrams of the original 
theory. The result is gauge fixing independent in covariant gauges and Coulomb-like 
gauges. Furthermore, it is transverse with respect to the external 4-momentum, which 
is necessary for the resulting effective theory to be gauge invariant. The dependence 
on the cutoff is logarithmic, and the cutoff dependent part was computed explicitly. 

In the physically interesting region of very small frequencies po -C g 2 T, only the 3-d 



17 



transverse part of the polarization tensor is unsuppressed relative to the hard thermal 
loop SU t . The longitudinal part is small compared to SUg , which, in the small frequency 
limit, approximately equals m^. 

There are higher loop diagrams in the hard thermal loop effective theory which are 
as important as the ones considered in this paper and which have to be included in the 
effective theory for the soft field modes. Evaluating them using the present framework 
is quite tedious. It is a lot more convenient |7|, |36|] to use the formulation of the hard 
thermal loop effective theory in terms of classical kinetic equations. In this approach 
it is relatively easy to resum all leading logarithmic contributions from higher loop 
diagrams. One obtains a Boltzmann equation for the field modes with p < fi. The 
diagrams computed in this paper together with all leading higher loop contributions are 
contained in a Gaussian white noise and a collision term in the Boltzmann equation. 
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Appendix A 

For diagrams containing the hard thermal loop resummed propagators, the sum over 



Matsubara frequencies is not easy to compute. The standard method |3T| is to write 
the propagators in their spectral representations. This method is not always useful. 
In particular, for the diagram in Fig. 2b it turns out to be very inconvenient. This 
Appendix describes an alternative method, which is very efficient when applied to the 
diagrams considered in Sect. 4. 
The sums are of the form 

S = T £ f(k Q ,p ), (Al) 

ko=iu> n 

where po = iuj n i is a Matsubara frequency. For given p , the function f(k ,p ) has cuts 
in the complex ko plane for Im^o) = and for Im(/c ) = — Po (Fig. ||). 

First we write the sum as an integral along the contour C a depicted in Fig. |5](a) 

S = J^(± + n(k ^f(k ,p ), (A.2) 

C a 



10 In scalar field theory the summation of ladder diagrams also leads to a Boltzmann equation 
containing a collision term J37|, 
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Im(fco) 



Re(fco) 











* 


3 






< 


) 







(a) 



(b) 



(c) 



Figure 5: Integration contours in the complex k plane for the integrals in (a) Eq. ( |A.2| ) 
and (b) Eq. ( |A.3|) . The contour in (c) corresponds to the integral in Eq. ( A. 6 ). The 



lowest and highest branches can be moved to ±ioo and do not contribute. The thick 
lines are the cuts of the function /(fco,Po) located at lm(k ) = and Im(A;o) = — lm.(p ). 



where n(ko) is the Bose distribution function (|3.1| ). Then the contour is deformed as 
in Fig. ||(b), which gives 

S = J^( K l + n(k )^(f(k + t5,p )-f(k -t5,p )) 

+ y Q + n(k ) \ (f(ko -po + iS, po) - f(k - Po - i6, Po)) ■ (A.3) 

Here the first (second) line corresponds to the upper (lower) part of the contour. In 
the second line we have used the periodicity of the function n, n(k +p ) = n(k ). Now 
we can perform the analytic continuation of po to arbitrary complex values. We are 
interested in the case fl2.9|). The k - integral is saturated for \k Q \ C T so that we can 
use the high temperature approximation 

\ + ni ~ ko) ~ V (A ' 4) 

to obtain 

S ^ J^jr(f(k + i5,p ) - f(k -i5,p )) 

+ J^jr(f{ko-Po + i8,Po)- f{ko-Pa-i8,p )), (A.5) 

where 5 < Im(po)- Now we deform the integration contour as in Fig. |5](c): The upper 
part is closed around the pole at k = 0. The piece above the cut can be moved towards 
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% x oo to give zero if f{ko,po) vanishes in this limit. The piece below the the cut is 
moved downwards to lm.(k ) = — Im(p /2). Proceeding similarly for the lower part of 
the contour, we finally obtain 

S ~ t(/(0,po) + /(-Po,Po) -Pa f^ T \ t \ f(ko ~ |pq,Po)1- (A.6) 
[ J 2m k + ±p k - ±p J 

Appendix B 

In this Appendix it is shown that one obtains the same leading order result for the 
polarization tensor as in Sect. |], if one uses Coulomb gauge rather than a covariant 
gauge. 

In Coulomb-like gauges the hard thermal loop resummed propagator reads 



K 2 h^h v 

w = A t vr + <fV° + H^r- 



fB.r 



Strict Coulomb gauge corresponds to £ — > 0. In order to see the effect of using the 
propagator ( |B.1| ) instead of ( |2.7p , it is convenient to write g^g 1 ' in terms of the tensors 
(cf. Ref. §) Vr, 



V2k 



9 



K 2 



K 2 



V 1 



and 



Using 



D 



k^k u 



fB.2) 



(B.3) 



ViD = 0, 



-Vi, 



(B.4) 



and 



tr(V e ) = tr(D) = tr(C 2 ) 
one finds 

£sT = A t vr + &e 



tr(C) = tr(V e C) = tr(CD) = 0, 



k k 2 



Wk 1 



k 4 



(B.5) 



(B.6) 



Now one can see that the propagator (|B.6| ) gives the same result as the propagator 
( f2.7| ). Each term in and D^ u contains at least one factor k^ or k u , which will be 
contracted with some velocity vector v. Due to the delta functions 5(v ■ K) all these 
terms vanish. 
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